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We study a type I see-saw scenario where the right-handed (RH) neutrinos, responsible for the 
light neutrino mass generation, lie at the electroweak scale. Under certain conditions, the strength 
of the charged (CC) and neutral current (NC) weak interactions of the Standard Model particles 
with the heavy RH neutrinos can be large enough to allow the production of the latter at the 
LHC, opening also the possibility of observing other low energy signatures of the new physics in 
the electroweak precision observables as well as in searches for rare leptonic decays or neutrinoless 
double beta decay. In this scenario the flavour structure of the indicated CC and NC couplings of 
the heavy RH neutrinos is essentially determined by the low energy neutrino parameters, leading 
to fairly strong correlations among the new phenomena. In particular, we show that the present 
bound on the /i — > e + 7 decay rate makes very difficult the observation of the heavy RH neutrinos 
at the LHC or the observation of deviations from the Standard Model predictions in the electroweak 
precision data. We also show that all present experimental constraints on this scenario still allow 
i) for an enhancement of the rate of neutrinoless double beta decay, which thus can be in the range 
of sensitivity of the GERDA experiment even when the light Majorana neutrinos possess a normal 
hierarchical mass spectrum, and ii) for the predicted fi — > e + 7 decay rate to be within the sensitivity 
range of the MEG experiment. 



I. INTRODUCTION 



One of the best motivated extensions of the Standard Model consists of introducing extra fermions that 
are singlets under the Standard Model gauge group, which we will call for dcfinitencss right-handed (RH) 
neutrinos. With this particle content, the Lagrangian contains a Yukawa coupling of the RH neutrinos 
with the left-handed leptonic doublets and the Higgs doublet, which leads to a Dirac neutrino mass when 
the electroweak symmetry is spontaneously broken. Besides, it contains a Majorana mass term for the 
RH neutrinos, which is a priori unrelated to the electroweak symmetry breaking scale. The most popular 
scenario in this model consists on assuming that the RH neutrino Majorana mass scale is much larger 
than the electroweak symmetry breaking scale, thus naturally leading to tiny Majorana neutrino masses 
through the renowned see-saw mechanism pQ. Furthermore, any other low energy effect of the RH neutrinos 
decouples at least with their mass squared, resulting in a tiny rate for the rare leptonic decays [2], a tiny 
leptonic electric dipole moment [2] and a tiny deviation of the electroweak observables from the Standard 
Model predictions [4-6J, in excellent agreement with presently existing experimental data. 

It should be borne in mind, however, that the mass of the RH neutrinos is a free parameter which can take 
any value between zero and the Planck scale. An interesting possibility arises when the RH neutrinos have 
masses in the range 0(1004- 1000) GcV. If this is the case, the new particles could be produced and detected 
at the Large Hadron Collider (LHC), if their couplings to the Standard Model particles are sizable 0|B]- 
This situation, at first sight bizarre in view of the tininess of the neutrino masses, can be realized in some 
well motivated scenarios, namely when lepton number is approximately conserved (see, e.g., [9lHljh More 
importantly, the contributions from the new particles to the low energy phenomena are no longer suppressed, 
opening the possibility of finding additional low energy signatures of the new physics with experiments at 
the intensity frontier. Furthermore, the existence of RH neutrinos with masses in the range of (100 -7- 1000) 
GeV with lepton number violating interactions can dramatically enhance the rate of neutrinoless double 
beta decay [TH [T3] , inducing rates that could possibly be observable at GERDA [T3] , even when the light 
neutrinos present a normal hierarchical mass spectrum. 

In this paper we analyze the constraints from various experiments on the scenario where the RH neutrinos 
are accessible to collider searches, as well as the interrelation between the different constraints. More 
specifically, we will derive the constraints that follow from the present bounds on the rate of the process 
fi — > e + 7, with a special emphasis on the relation to the neutrino mixing and oscillation parameters. We 
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will also discuss the prospects to observe neutrinoless double beta decay in the next round of experiments, 
in the view of all present experimental constraints. In section 2 we discuss the formalism and define the 
parameter space of the theory. In sections 3 and 4 we discuss the constraints on the parameter space which 
arise from radiative charged lepton decays and the implications for collider searches and electroweak precision 
observables. In the subsequent section we perform a detailed analysis of the possible enhancement of the 
neutrinoless double beta decay rate due to the exchange of the heavy (RH) Majorana neutrinos. In section 
6 we combine the constraints on the parameter space that we obtained in sections 3, 4 and 5. Finally, we 
report our conclusions. 

II. PRELIMINARY REMARKS 

We denote the light and heavy Majorana neutrinos with definite masses as Xi an d N k , respectively. 1 The 
charged and neutral current weak interactions involving the light Majorana neutrinos have the form: 

Ccc = -4= h a "a. W a + h.c. = -A= ha ((1 + V)U) U XiL W a + h.c. , (2.1) 

Cnc = -7^WLl a vi L Z a =-^x^l a (uHl + V + rf)U) x 3 LZ a , (2.2) 

where (1 + rj)U = C/pmns is the Pontecorvo, Maki, Nakagawa, Sakata (PMNS) neutrino mixing matrix 
[ToT - fTH] . U is a 3 x 3 unitary matrix which diagonalises the Majorana mass matrix of the left-handed (LH) 
flavour neutrinos vu, (generated by the see-saw mechanism), and the matrix rj characterises the deviations 
from unitarity of the PMNS matrix. The elements of C/pmns are determined in experiments studying the 
oscillations of the flavour neutrinos and antineutrinos, vg and Dg, £ — e, t, at relatively low energies. In 
these experiments the initial states of the flavour neutrinos, produced in some weak process, are coherent 
superpositions of the states of the light massive Majorana neutrino Xi only. The states of the heavy Majorana 
neutrino Nj are not present in the superpositions representing the initial flavour neutrino states and this 
leads to deviations from unitarity of the PMNS matrix. 

The charged current and the neutral current weak interaction couplings of the heavy Majorana neutrinos 
Nj to the W and Z° bosons read: 

Lie = - ^= ha (RV)i k (l - 75) N k W a + h.c. , (2.3) 
C% c = -^-WZl a {RV)g k (l- 1 r )N k Z a + h.c. (2.4) 

Here V is the unitary matrix which diagonalises the Majorana mass matrix of the heavy RH neutrinos and 
the matrix R is determined by (see [T3]) R* — Md , Mjj and Mjv being the neutrino Dirac and the 
RH neutrino Majorana mass matrices, respectively, \Md\ <§C The matrix 77 which parametrises the 

deviations from unitarity of the neutrino mixing matrix, can be expressed in terms of the matrix R: 

V = ~RI$ = -\(RV)(RVy = r?t . (2.5) 

It is possible to constrain the elements of the hermitian matrix rj by using the existing neutrino oscillation 
data and data on electroweak (EW) processes [SJ E] (e.g., on W ± decays, invisible Z decays, universality 
tests of EW interactions). For M k k5 100 GeV these constraints read 0|6]: 

/ 2.0 x 10~ 3 0.6 x 10~ 4 1.6 x 10~ 3 \ 
\r)\ < 0.6 x 10~ 4 0.8 x IO- 3 1.0 x 10" 3 . (2.6) 

y 1.6 x io~ 3 l.o x io~ 3 2.6 x io~ 3 J 

The elements of the matrix RV and the masses M k of the heavy Majorana neutrinos Nk should satisfy 
the approximate constraint on the elements of the Majorana mass matrix of the LH flavour neutrinos [19 , 



1 We use the same notations as in I15| . 



3 



Km,,)^! < 1 eV, £, £' = e, n, t. In the case of the type I see-saw mechanism under discussion this implies: 

J2\(RVr e , k M k (RV)l\<leV, £',l = e,n,r. (2.7) 

This relation can be satisfied in several situations. The most trivial way to satisfy it is by imposing that 
| (-RV)^fe | <€. 1 for all £ and fc, which renders the observation of all RH neutrinos impossible at the LHC or 
in any low energy phenomena. However, this relation can also be fulfilled if one clement of the matrix RV 
is sizable. This requires the existence of at least another large matrix element, in order to cancel the large 



contribution of the former one in eq. (2.7 1. Whereas the possibility of cancellations between three different 



terms cannot be precluded, in the simplest case only two terms will cancel. Then, the large contribution to 



eq. (2.7 1 from one of the RH neutrinos, say Ni with mass Mi, is cancelled by a negative contribution from 



another RH neutrino, say N 2 with mass M 2 , provided 

[Aft 

(RV)e2 = ±i(RV)a\ -± (2.8) 

V J- VI 2 

which is naturally fulfilled if the RH neutrinos N\ and N 2 form a pseudo-Dirac pair [501 [21]) i-e., if there 
exists an approximately conserved lepton charge (see, e.g., [H]). In this scenario, in order not to spoil 



the cancellation between these two terms, the contribution from the third neutrino to eq. (2.7) should be 
negligible. Therefore in what follows we will work for simplicity in the 3x2 see-saw scenario, in which the 
indicated CC and NC weak interaction couplings of N$ are set to zero and N$ is decoupled. In this case the 
SM is effectively extended by the addition of 2 RH neutrino fields only. In this class of models (see, e.g., 
[22Tf24] ) one of the three light (Majorana) neutrinos is massless and hence the neutrino mass spectrum is 



hierarchical. Two possible types of hierarchical spectrum are allowed by the current neutrino data (see, e.g., 
[2"5]): i) normal hierarchical (NH), mi = 0, m 2 = J Am| and 7773 = y'Am} , where Atoq = m\ — m 2 > 
and Am\ = m| — to 2 ; 



ii) inverted hierarchical (IH), 7773 = 0, 771,2 = \J\Am\ I an( i m i = \/l^ m A I — ^ m © — \/l^ TO A I' wnere 



Atoq = ml — m\ > and Am\ = m\ — m\ < 0. In both cases we have: Atoq/IAto^ | = 0.03 <C 1. 

T 

The two heavy neutrino fields satisfy the Majorana condition: CN k = N k , fc = 1,2. If the Majorana 
mass matrix M/v of the RH neutrinos is not CP invariant, one can always make the eigenvalues of Mn real 
and positive, M\$ > 0. For Mjy respecting the CP symmetry, the two real eigenvalues of Mm can have 
the same or opposite signs (see, e.g., [26]). With the choice M 2 > we can always make (without loss of 
generality), Mi can be positive or negative: M\ > and Mi < 0. One can show, however, that we get the 
same results for the observables of interest to this study in the two cases. Therefore in what follows we shall 
work with Mi. 2 > 0. 

We assume that the heavy RH neutrinos 7V li2 have masses in the range Mi 2 = 0(100-7- 1000) GeV, which 
makes possible, in principle, their production, e.g., at LHC. In order to be produced with observable rates 



at LHC, the CC and NC couplings of N% and N2 in eqs. (2.3) and (2.4) have to be sufficiently large. Under 
this condition the existing experimental upper bounds on the neutrinoless double beta ({fift)ov) decay rate 
put stringent constraints on the mass spectrum of the RH neutrinos. It can be shown [T5], in particular, 
that in the type I see-saw scenario of interest the two heavy (RH) Majorana neutrinos must be almost 
degenerate in mass: M 2 = Mi. If we assume that M 2 > Mi > and M2 = (1 + z)M\, z > 0, in order 
to satisfy the experimental limit on the (/3/3)oi/-decay rate, one should have z < 10~ 3 (10~ 2 ) for Mi ps 10 2 
(10 3 ) GeV [15] . 

The charged current and the neutral current weak interaction couplings of the heavy Majorana neutrinos 
Nj, (RV)gk, are furthermore constrained by the requirement of reproducing the correct low energy neutrino 
oscillation parameters after the decoupling of the heavy degrees of freedom. Remarkably, under the condition 
that \(RV)ik\ are sufficiently large to produce observable effects of the RH neutrinos at low energies, these 
couplings take a very simple form [15] . 

Indeed, the Dirac mass matrix Mq in the case under study can be written as |27j 

Md = jL^ MJVS Vm"oVMF t , (2.9) 
where 777 = diag (to 1, 7772, 7773) and O is a complex orthogonal matrix. In the scheme with two heavy RH 
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Majorana neutrinos the matrix O has the form 



O 



O 





cos 9 ± sin 9 
— sin 9 ± cos 9 



sin ( 




± sin 9 
± cos 9 




for NH mass spectrum 



for IH mass spectrum 



(2.10) 
(2.11) 



where 9 = u — i^. The RH neutrino mixing matrix entering into the CC and NC weak interaction Lagrangians 
(2.3) and (2.4) can be expressed as [15] 



RV = -i Upmns VmO' V M" 1 . 
The O-matrix in the case of, e.g., NH spectrum can be decomposed as follows: 

o = i 





(2.12) 



(2.13) 



One can get a similar expression for the IH spectrum. The Dirac neutrino mass matrix can be decomposed 
accordingly as = Mo+ + Md-i in a self-explanatory notation. Taking for definitcness £ > 0, it follows 
that Mo + {Mo-) grows (decreases) expon entia lly with £. 2 Therefore, for sufficiently large £ it is possible 
to compensate the huge suppression in eq. (2.9) from the tiny observed neutrino masses and the relatively 
light RH neutrino masses. 3 

We are interested in heavy Majorana neutrino couplings to charged leptons and gauge bosons, which are 
large enough to produce observable low energy signatures. In the limit of "large" £, the matrix O in eq. 
(2.13) can be very well approximated by: 



O 



For NH spectrum and "large" £, the matrix RV takes the form: 




(2.14) 



RV 



m 3 
|M X | 



^ (U e 3 +iy/m 2 / m 3 Ue2) ±i (u e3 + i^/m 2 /m 3 U e2 ^j j\f\ + z ^ 
y (U T 3 + iy/m 2 /m 3 U r2 



±i IU ( 

±i + i\Zm 2 /m 3 / y/T+ 
±i (Ut3 + i\/m 2 /m 3 U r2 j + Z , 



(2.15) 



In a different context similar expressions for the corresponding neutrino Yukawa couplings were derived 
in [28] in a scheme in which a successful leptogenesis can take place at T ~ 10 7 GeV, and in [11] where 
a TeV scale see-saw model with approximately conserved lepton charge was proposed. In the case of IH 
light neutrino mass spectrum, the matrix RV is obtained by replacing m 2 3 —¥ m 12 and U a2 ^ a 3 — > U a \^ a2 
(a = e,/i,r) in eq. (2.15). For both types of neutrino mass spectrum, the elements of the two columns of 



RV in the limit of "large" £ are related by the condition: 

(RV) a2 = ±i (RV) al /Vl + z for a = e,fi, r. 



(2.16) 



2 Obviously, if £ < 0, Mo— (Mn+) will grow (decrease) exponentially with £. It is possible to show that for sufficiently large 
values of |£| of interest, the results for the different observables considered in our study do not depend on the choice of the 
sign of £. 

3 Note, however, that Mq_ cannot be neglected in the calculation of the Majorana mass matrix of the LH flavour neutrinos 
even though it is exponentially suppressed compared to Mrj>+: the naive approximation Mjj ~ leads to m v = 0, due 
to O+O^ = (see |15| for details). Therefore, reproducing the correct light neutrino masses and mixing requires a large 
amount of fine tuning, unless the RH neutrinos form a pseudo-Dirac pair. Demanding [Mo)ij ~ 0(1 GeV) and random 
RH neutrino masses of 0(100 GeV), for instance, requires a tuning of one part in 10 9 in order to produce a neutrino mass 
in, ~ O(10" 2 eV). 
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thus recovering eq. (2 



If Mjv is CP conserving an d Mj < 0, Ma > 0, the expressions for [RV) a \ in eq. (2.151 have an additional 
factor —% an d ins tead of eq. (2.16) we have: {RV) a 2 — T (W) a i/Vl + z - This relation and the relation 
given in eq. ( |2.16 1 lead to the same expressions for the observables discussed further in the present study. 

The overall size of the couplings eq. (2.15 1 depends crucially on the value of the parameter £ which has no 



direct physical interpretation. It proves convenient to express £ in terms of the largest eigenvalue y of the 
matrix of neutrino Yukawa couplings using the relation: 



2 2 

V v 



jeig (McA/jj) } = max jeig^N/^OMC^V™)} = -e 2? Mi(m 2 + m 3 )(2 + z) , (2.17) 



with v = 174 GeV. In terms of y and for z <C 1, the heavy Majorana neutrino couplings become: 



I Wall 
I Wall 



1 y 2 v 2 



m 3 



2 M{ m 2 + m 3 

7712 



1 y 2 v 2 



2 M 2 mi + m 2 



U a z + i\/m 2 /m 3 U a2 
U a2 + iyj mi/m 2 U a i 



, NH, 

2 1 „,2„,2 



1 V V > 

~r\u a2 + iu al \ 2 ,ih, 



(2.18) 
(2.19) 



where we have used the fact that for the IH spectrum one has mi = m2- 

From eqs. (2.18) or (2.19), using the unitarity of the matrix U, one can express the neutrino Yukawa 



eigenvalue y in (2.17) in terms of heavy neutrino to charged leptons coupling constants: 

fl (iWel 



2 2 
V V 



2 Mi 



IW 



/' i 



IWrl 



(2.20) 



An upper limit on the neutrino Yukawa coupling y can be derived by assuming the validity of perturbative 
unitarity. Indeed, the requirement of perturbative unitarity can be easily fulfilled if the see-saw parameter 
space satisfies the condition: 



M, 2 



(2.21) 



where is the total decay rate of the heavy Majorana neutrino N t , which in the limit Mjf ( 3> v is given 
by the following expression: 
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WttM 2 , 



Mf £|(W)«| 2 , 



(2.22) 



limit on the parameter y from perturbative unitarity: 

y < 4 



Therefore, taking into account eqs. ( |2.20[ ) and ( |2.22[ ) and the condition (2.21), we get the following upper 

(2.23) 



III. THE fj, -»• e + 7 DECAY 

The fi — > e + 7 decay branching ratio in the scenario under discussion is given by [291 130] : 



B{n ->• e + 7) 



T{n e + 7) 



3 a 



— >• e + z/ M + i/ e ) 327T 



cm \rp\2 



where a rm is the fine structure constant and 



(3.1) 



3 = 1 
2 



( M 2 



k=l 



(3.2) 



() 



The loop integration function G{x) has the form: 

10 - 43a; + 78a; 2 - 49x 3 + 4x 4 + 18a; 3 log(x) 



G(x) 



3(x-l) 



(3.3) 



It is easy to verify that G(x) is a monotonic function which takes valu es in the interval [4/3, 10/3], with 

10 
3 



G{x) = ^ — x for x <C 1. From the definition of the matrix r], eq. (2.5), one can write the amplitude T as 



follows: 



T = [(RV);, (RV) el + (RV); 2 [RV) e2 ] [G(X) - G(0)] 



(3.4) 



where X = (M\ / Mw) 2 and we have assumed that the difference between Mi and Mi is negligibly small and 
used Mi = M 2 - Using eqs. (3.4) and (2.16) we get: 



in s 



1 + 2 



\(RV)* (RV) el \\G(X)-G(0)\ . 



(3.5) 



Finally, using the expressions of K-RU)^! 2 and |(-RV%i| 2 in terms of neutrino parameters, eqs. (2.18) and 
(2.19), we obtain the fi — > e + 7 decay branching ratio for the NH and IH spectra: 

B(fi ->• e + 7) = 



NH : 



2 2 
7/ z 7T 



m 3 



32-7T \ Mf m 2 + m 3 



+ H\'—Uva 
u m 3 ' 



t/ e 3 + «4 



-0; 



c2 



[G(X)-G(0)] s 



IH: fl(/j->e + 7) 



32tt V M- 



2 2 1 \ 2 
j/ z tr l 



T^) 1^2 + «^i| 2 |C/ e2 + iC/eil 2 [GpO - G(0)] 



(3.6) 



(3.7) 



Employing the standard parametrisation of the neutrino mixing matrix |25j it is not difficult to obtain 
expressions for the factors |[/^ 3 + i ^/ '777,2 / 'm^U £ 2 | 2 and \Un + iUn\ 2 , £ = e, fj,, in terms of the neutrino mixing 
parameters and the solar and atmospheric neutrino mass squared differences: 




7772 

U e .3 + i\j — U e2 

777 3 



'13 



Am 2 n 1/2 



Am\ 



c 13 s 12 



A777q 

Ami 



1/4 



ci 3 si 3 S12 sin 5 



«2i - « 3 i 



(3.8) 



C^3 + Hl'—Unz 

W 777 3 



c 13 s 23 



1/2 



( c 12 c 23 + s 12 s 23 s 13 ~ 2c 12 C23 Sl2 S23«13 COS 5) 



+ 2 



/ Am|, 



Cl3 S 2 3 



Ci 2 C 23 sin 



c* 3 i - a 2 i 



Sl2 S 23 Sl3 Sill 



C*31 - Ct 2 l 



13 



1 + 2c i2 S12 sin {^y) 



\Uft2 + iU^l | = C 2 3 + S 2 3 S2 3 - 2 C12 S12 {c\ 3 - Si 3 s\ 3 ) SU1 ■ 



- 2 C 23 Si 3 S23 



. / "21 , A 2 . /«21 A 

: Bin ^— + <SJ - c 12 sin ^— - <5 j 



(3.9) 
(3.10) 

(3.11) 



where c, 



cos 9ij, Sij = sin 9ij, 12 , #23 and #i 3 are respectively the solar neutrino, atmospheric neutrino 



and the CHOOZ angles, S is the Dirac CP violating phase and a 2 i and a 3 i are the two Majorana CP 
violating phases (3"T] . 

We show in Fig. [I] the branching ratio of \i — > e + 7 as a function of the RH neutrino mass Mi , for three 
different values of the neutrino Yukawa coupling eigenvalue: y = 0.001 (0.01) [0.1], blue o (green +) [red x]. 
The figure was obtained by performing a scan of the values of the neutrino mixing angles #i 2 , # 23 and 813 
and the solar and atmospheric neutrino mass squared differences At77q and Am\ within the corresponding 
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3<7 bounds (see Table [T]). 4 The Majorana phases a^i and (0:31 — a.21) are varied in the interval 5 [0, 47r] [33] 
and the Dirac phase 5 is varied in the interval [0, 2tt]. Shown are also the current experimental upper limit 
on the /i — > e + 7 decay branching ratio [33], B(/i — > e + 7) < 1.2 x 10 -11 , as well as the prospective limit of 
the MEG experiment [33], B(/i ->• e + 7) < 10" 13 . 

It is apparent from Fig. [l] that the data on the process [i — > e + 7 set very stringent constraints on the 
TeV scale see-saw mechanism. A neutrino Yukawa coupling y = 0.1 generates a rate for the \x — > e + 7 decay, 
which is ruled out by the MEGA experiment, unless the RH neutrino mass is M\ > 300 GeV. Furthermore, if 
the MEG experiment reaches the sensitivity of 10 -13 without finding a signal, for the same Yukawa coupling 
y — 0.1 the RH neutrino mass should be larger than 1 TeV. More generally, for Mi = 100 GeV (M\ — 1 
TeV) and z < 1 we get the following upper limit on the product \(RV)* fll (RV) e i\ of the heavy Majorana 
neutrino couplings to the muon (electron) and the W boson and to the Z° boson from the current upper 
limit [33] on B(/i — > e + 7): 



\{RV)* Ml (RV) e i\ < 1.8 x 10~ 4 (0.6 x 10~ 4 ) , 



(3.12) 



where we have used eqs. (3.1) and (3.5). This can be recast as an upper bound on the neutrino Yukawa 
coupling y. Taking, e.g. the best fit values of the solar and atmospheric oscillation parameters, we get: 



y < 0.036 (0.21) for NH with M x = 100 GeV (1000 GeV) and sin 13 = , 
y < 0.031 (0.18) for IH with M x = 100 GeV (1000 GeV) and sin0i 3 = 0, 
y < 0.094 (0.54) for NH with M 1 = 100 GeV (1000 GeV) and sin0i 3 = 0.2 . 
y < 0.16 (0.90) for IH with M x = 100 GeV (1000 GeV) and sin6> 13 = 0.2. 



(3.13) 
(3.14) 
(3.15) 
(3.16) 



The upper limit on the neutrino Yukawa coupling derived for #13 = is reached for 021 — a.31 — 7r (0:21 — 37r) 
in the case of NH (IH) light neutrino mass spectrum. For sin 13 = 0.2, the upper limit on y is obtained for 
c*2i — CK31 — I" (&21 — t) an d 8 ~ (S ~ 0), if the neutrino masses have normal (inverted) hierarchy. 

A possible way to circumvent these stringe nt up p er bou nds consists in assuming a very fine cancellation 
between the different terms in one of th e eqs. (|3 .8 ) — (3.11 ) for the factors \Ut3 + zy 'milm^Un^ and \Ui<x + 

7 ) for B(n — > e + 7) . Such a cancellation is possible in the 
e3 + i^/ r m 2 /m 3 U e2 \ = if 



iU ei \ 2 , I = e,n 

case of NH spectrum (see also 



in the expressions (3.6) and 

^8]). Indeed, we have \U, 



sm 



Q?21 ~ Q 31 

2 



= 1, and tan 0i3 = 



Am% 



1/4 



sin P12 



(3.17) 



In this case B(p — > e + 7) = and, thus the upper bound on the Yukawa coupling y is no longer valid. 
Assuming the neutrino parameters have values within the corresponding present 3cr ranges reported in 
Table 1, we find that the values of 0i 3 implied by the condition in eq. ( 3.17| satisfy: tan0i3 > 0.21, or 
sin 2 0i3 > 0.043. These values are outside the 3tr interval of the experimentally allowed values of sin 2 0i3 
(see Table 1). 

We consider next the case in which the neutrino mass spectrum is with inverted hierarchy. The first thing 
to notice is that the factor \U e 2 +iU e i\ 2 can be rather small for sin(a2i/2) = —1 since 2ci2 S12 — 0.93, where 
we have used the best fit value of sin 2 i2 = 0.312. In this case we have \U e 2 + iU e \\ 2 = 0.069 (0.066) for 
sin 2 0i3 = (0.04). Second, we show that, as in the case of normal hierarchical spectrum, it is possible to 
have a strong suppression of the factor \U^2- s riU ll i\ 2 if the CHOOZ mixing angle is close to the corresponding 
3cr experim ental upper bound. To be more quantitative, we take sin 2 0i2 = 1/3 and sin 2 023 = 1/2. Then, 
expression (3.11) takes the form: 



I[ 2 (V2(-1 



s 13) ~ s 13 COS S 



sm ■ 



«21 



3 1 



'13 



„ "21 . 

2si 3 cos — - sm 



(3.18) 



4 In all scatter plots included in our paper the neutrino observables are scanned assuming a gaussian distribution with the 
corresponding mean value and standard deviation reported in Table 1. All the other (unmeasurable) see-saw parameters are 
selected with a flat distribution. 

5 We note that the phases «2i and «3i enter into the expression for the neutrino mixing matrix in the form exp(io2l/2) and 
exp(ia3i/2), respectively. 



Normal Hierarchy 



Inverted Hierarchy 



ir 10- 4 



i, 10 




100 150 200 300 500 700 1000 15002000 

M\ (GeV) 



100 150 200 300 500 700 1000 15002000 

M\ (GeV) 



FIG. 1: The dependence of B(fi — > e + 7) on Mi in the case of NH (left panel) and IH (right panel) light neutrino 
mass spectrum, for i) y = 0.001 (blue o), ii) y = 0.01 (green +), and Hi) y = 0.1 (red x). The horizontal dashed line 
corresponds to the MEGA bound [S3], B(fi — > e + 7) < 1.2 x 10 -11 . The horizontal dot-dashed line corresponds to 
B(/i — > e + 7) = 10 -13 , which is the prospective sensitivity of the MEG experiment |34| . 



It is not difficult to show that, for fixed values of the phases «2i and 8, |I7 M 2 + *t^ii| 2 has a minimum for 



sin 6*13 



cos 8 sin 



3 cos ^ sin 8 



3 + 2^2 sin- 



(3.19) 



At the minimum, using eqs. (3.181 and (3.19), we get: 



M\ u > 



//2 



iU„i\ 2 ) 



(3 cos 8 cos ^ + sin 8 sin 
6(3 + 2\/2sm^) 



(3.20) 



We will find next for which values of the CP violating phases 8 and a.21 this lower bound is equal to zero 
and if the resulting # 13 , obtained from eq. (3.19), is compatible with the existing limits from the neutrino 
oscillation data. We have min(|C/ At 2 



if the Dirac and Majorana phases 5 and 021 satisfy 



the following conditions: tan^tan^ 1 = —3 and sgn(cos 8 cos 9 ^ L ) = — sgn(sin<5sin 9 ^ L ). Taking cos<5 > 



(cos 5 < 0) and using tan<5 = —3/ tan(a2i/2) in eq. (3.19) we get: 



sin 6*i3 = sgn(cos<5) 



9 + tan 2 S|i 



3 + 2\/2sin- 



C*21 



(3.21) 



The solution (3.21 



) is compatible with the 3cr upper limit of the CHOOZ mixing angle (see Table 1). In 
general, one can always find a viable pair of CP violating phases a2i and 8 satisfying the relations given 
above in order to set the r.h.s. of eq. (3.20) equal to zero, if the mixing angle #13 is sufficiently large, namely, 
if sin 6»i3 > 3 - 2^/2 = 0.17. More precisely, one finds, e.g. that \U^ 2 + iU^\ 2 ~ 3.52 x 10~ 8 (2.43 x 10~ 6 ) 
for S13 ~ 0.2 (0.17), a 2 i ~ 2.732 (tt) and 8 ~ 5.725 (10~ 3 ). 

In order to interpret the results presented in Fig. [T] it proves convenient to use the analytic expressions 
of B(p — > e + 7) in term s of th e low energy neutrino parameters, the neutrino Yukawa coupling and the 
RH neutrino mass, eqs. (3.6 1— (3.11 ). Taking for concreteness sin 2 #23 = 1/2, sin 2 6*12 = 1/3 and using 
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sin 0i3, Am s 2 |/Am3 tm <§C 1, eqs. (3.6) and (3.7) approximately read 



NH : B(fi -> e + 7) 



-> / 2 2 \ ^ 
3a nm / y z w z 



32vr 



Mi 2 



[G(X)-G(0)] s 



6 V Am\ 



1/4 



(3.22) 



/ Ami 
VAmT 



IH : #(/z 



1/4 



2\/3si 



sin 0i3 sin (5 



«21 - «3i 



2 / Am| 



1/4 



0:21 - «31 



, 3a 
■7) - " 



32tt 



y 2 v 2 1 
M^2 



[Gpf) - G(0)] 2 x 



1 
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5 + 4cosa2i — 2 sin 0i3 (3 + 2y2sin 



a2i 



)( 



cos 5 sin 



3 sin S cos 



«21 



(3.23) 



7) 



From these expressions it follows that there is a fairly strong dependence of the prediction of B(/j, — > 
on the Majorana phases. As a result, and since there is no proposal to constrain experimentally these phases, 
there is uncertainty band of a factor of five for normal hierarchy and a factor of nine for inverted hierarchy, 
which cannot be reduced even if the neutrino masses and mixing angles were known with arbitrarily high 
precision. 

It is not difficult to get also expressions for the double ratios #(21/31) = #(/x — > e + j)/B'(t — > e + 7) 
and #(21/32) = B{ji e + j)/B'(t (J, + 7), where B'(t -> e(/i) + 7) = B(t -> e(/i) + j)/B(t ->■ 
e(/i) + v T + v e (^))i B(t — > e(p) + 7) and B(j — » e(/x) + ^ T + S'eO)) being the branching ratios of the 
corresponding decays. The double ratios of interest depend only on the neutrino masses and the elements 
of the PMNS matrix: 



#(21/31) 
#(21/31) 



^3 + H/^2 



m 3 
/ m 2 , 



^3 + ^/^2 



NH. 



I^r2+it/rl| 2 



IH: 



(3.24) 
(3.25) 



#(21/32) 
#(21/32) 



U e3 + i 



V2 



TO -W m3 TZ 



NH. 



IK 



e2 



■ill, 



el 



\U T2 +iU Tl [ 



, IH. 



(3.26) 
(3.27) 



For the NH (IH) light neutrino mass spectrum, the range of variability at 3<r of each of the two ratios defined 
above is: 0.01 < #(21/31) < 20 (0.001 < #(21/31) < 300) and 5 x 10~ 4 < #(21/32) < 3 (0.04 < #(21/32) < 
5000). 6 Thus, in the case of NH spectrum, the predicted r — > \i + 7 decay branching ratio B'(r — > /1 + 7) 
in the scheme considered can be by several orders of magnitude larger than the fi — > e + 7 decay branching 
ratio, B(fi — > e + 7). 

IV. IMPLICATIONS FOR COLLIDER SEARCHES AND ELECTROWEAK PRECISION 

OBSERVABLES 

Upper bounds on the couplings of RH neutrinos with SM particles can be set by analyzing lepton number 
conserving processes like ir — > £~p£, Z —¥ vv and other tree- level processes involving light neutrinos in the 



In order to get such estimates we assume that the neutrino observables have a gaussian distribution with the corresponding 
mean value and standard deviation reported in Table 1. 
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Parameter 


riest r it ± la 


ocr interval 




NH | IH 


NH | IH 


Amlj (HT 5 eV 2 ) 


7 50+0-20 


[7.09,8.19] 


lAmlxl (10" 3 eV 2 ) 


9 ^r+0.09 9 o 4 +0.10 
z -^ J -0.09 z -°^-0.09 


[2.18,2.73] [2.08,2.64] 


sin 2 612 


n oi 9+0.017 


[0.27,0.36] 


sin 2 823 


0.51 ±0.06 0.52 ±0.06 


[0.39,0.64] 


sin 2 6 13 


1 0+ 009 nn+ 009 
u.uiu_ 006 u.uio_ 007 


< 0.035 < 0.039 



TABLE I: Best fit values with la errors and 3a intervals for the three flavour neutrino oscillation parameters (see 
[35| and references therein). 



i-3 



final state J5]- From eqs. (2.5) and (2.6), we get: 

|(W) el | 2 < 2x 10" 

IW^il 2 £ 0.8x io- 3 : 

l(-R^) r i| 2 < 2.6 x 10~ 3 . 



(4.1) 
(4.2) 
(4.3) 



Following the same rationale as in the previous section, one can translate the upper bounds on the RH neu- 
trino couplings from electroweak precision observables into upper bounds on the neutrino Yukawa coupling 
y. For this purpose we observe that an upper limit on y, which is independent on the specific values of the 
neutrino oscillation parameters, can be easily derived from eq. (2.20), taking into account the experimental 
constraints given above: 



V < 0.06 



( Mx 



V 100 GcV 



(4.4) 



On the other hand, using eq. (2.19) and taking the best fit values of neutrino oscillation parameters, we 
obtain 



y < 0.047 (0.47) for NH with Mi 
y < 0.046 (0.46) for IH with M x = 
y < 0.049 (0.49) for NH with M x 
y < 0.053 (0.53) for IH with M x = 



= 100 GeV (1000 GeV) and sin 13 = , 
100 GeV (1000 GeV) and sin U = , 

= 100 GeV (1000 GeV) and sin6>i 3 = 0.2, 
100 GeV (1000 GeV) and sin 13 = 0.2 . 



(4.5) 
(4.6) 
(4.7) 
(4.8) 



The results reported in eqs. ( |4.5[ ) and ( |4.6| ) arc obtained for a 2 \ — «3i — 7r and a 2 i — 0, respectively. In 
the case of sin #13 = 0.2 and NH (IH) light neutrino mass spectrum, the upper bound on y corresponds to 
c*2i ~ a 3i — n i a 2i — 0) and 5 ~ ir (S ~ 37r/2). 

The bounds for sin ^13 = are weaker than those derived in the previous section from the non-observation 
of the fx e + 7 decay. Thus, the existing stringent upper bound on the /x — >• e + 7 decay rate makes very 
difficult the observation of deviations from the Standard Model in the electroweak precision data, predicted 
by the TeV scale see-saw scenario. In contrast, in the case of sin 6*13 = 0.2, the constraint from the MEGA 
upper bound [33] can be avoided and we get a better upper limit on y from the electroweak precision data. 
More precisely, we find that for Mi = 100 (1000) GeV, the limits given in eqs. ( |4. 1| ) — ( [473] ) provide a better 
constraint on the neutrino Yukawa coupling y than the upper bound on B(/i — > e + 7) if sin 6*13 > 0.10 (0.19) 
in the case of the NH neutrino mass spectrum, and for sin #13 > 0.13 (0.17) if the spectrum is of the IH type. 
In particular, the parameter y can have a value as large as 47r that is still compatible with the current upper 
limit on B(fi e+7) [33]. This is possible in the case of NH (IH) light neutrino mass spectrum for sufficiently 
large values of sin6»i 3 > 0.22 (0.18). We note, however, that at sinfli 3 = 0.10 (0.19) (sin0i 3 = 0.13 (0.17)) 
for the NH (IH) spectrum and M x w 100 (1000) GeV, we have y ^ 0.05 (0.5) both from the bound on 
B(ji — > e + 7) and the electroweak data limits. This is consistent with the absolute upper limit reported in 
(ph. 
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FIG. 2: Correlation between |(-RV) e i| and |(J?V) M i| in the case of NH (left panels) and IH (right panels) light 
neutrino mass spectrum, for Mi = 100 (1000) GeV, upper (lower) panels, i) y = 0.001 (blue o), ii) y = 0.01 (green 
+), Hi) y = 0.1 (red x) and iv) y = 1 (orange (>)■ The cyan points correspond to random values of y < 1. The 
dashed line corresponds to the MEGA bound [33], B(fi — > e + 7) < 1.2 x 10 -11 . The dot-dashed line corresponds to 
B(/i — > e + 7) = 10~ 13 , which is the prospective sensitivity of the MEG experiment |34| . 



These results are illustrated in Fig. [2j where we show, for Mi = 100 GeV (upper panels) and Mi = 1000 
GeV (lower panels), the allowed ranges of the RH neutrino couplings |(i?V) M i| and |(i?V) e i|, in the case of 
NH (left panels) and IH (right panels) spectra. The region of the parameter space which is allowed by the 
electroweak precision data, eqs. (4.1 ) and (4.2 ), is indicated with solid lines; the region allowed by the current 
bound on the /1 — s- e + 7 decay rate is indicated with a dashed line. The projected MEG sensitivity reach 
is shown with a dot-dashed line. Finally, in the same figure we show a scatter plot of the points which are 
consistent with the 3a allowed ranges of the neutrino oscillation parameters. This is done for four different 
values of the neutrino Yukawa coupling y: y = 0.001 (blue o), ii) y = 0.01 (green +), in) y = 0.1 (red x) 
and iv) y — 1 (orange ()). The cyan points in each plot correspond to an arbitrary value of the Yukawa 
coupling y < 1. 

After imposing all the constraints on the parameter space discussed above, we find a fairly narrow band 
of values of |(i?V) M i| and |(i?V) e i|, which is allowed by the data, more precisely, by the requirement of 
reproducing the correct values of the neutrino oscillation parameters and by the constraint following from the 
upper bound on the ji — > e + 7 decay rate. Interestingly, and as we have already mentioned earlier, the limits 
from the electroweak precision data lie essentially in the excluded (shaded) region. Therefore, in the view of 
the constraints from the data on the neutrino oscillation parameters and the /1 — > e + 7 decay, the discovery of 
significant deviations from the Standard Model predictions in the electroweak precision observables appears 
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highly improbable, unless a significant improvement in the precision of the data is achieved. The above 
conclusion will be strengthen if the MEG experiment reaches the sensitivity B(fi — > &y) ~ 10~ 13 without 
observing a signal. 

It has been argued that present constraints from electroweak precision data allow the observation of 
pseudo-Dirac RH neutrinos at the LHC. More concretely, the following process with three charged lepton 
final state, 



qq -> n r N PD -> W + -> ^ + \i /x+z^ . (4.9) 

could be observed at the LHC with a luminosity of 13 fb _1 with a 5<r significance if the RH neutrinos have 
a mass of ~ 100 GeV [5]. A discovery reach in this channel implies a coupling |(iJV")ui| ~ 0.04 [5]. Taking 
\ £ 0.04, we get from ( pUf and ([2T9|: 

= 100 GeV, (4.10) 
100 GeV. (4.11) 



y > 0.04 for NH with Mi 
y > 0.05 for IH with M x = 



which should be compared with the upper limits on y we get from the upper bound on the /i — > e + 7 decay 
rate, eqs. (3.13) and (3.14), or from electroweak precision data, eqs. (4.7) and (4.8). We find again a tension 
between the constraints on y obtained from the data on the /1 — > e + 7 decay or the electroweak processes 
and the values of y required for the production of RH neutrinos with observable rates at colliders. This is 
clearly seen in Fig. [2j 

As a conclusion, the presently existing data on the neutrino mixing parameters and the present experi- 
mental upper bound B(fi —> e + 7) < 1.2 x 10~ n , basically rule out the possibility of producing pseudo-Dirac 
neutrinos at LHC with observable rates. A similar conclusion applies to the possibility of observing devia- 
tions to the Standard Model predictions in the electroweak precision observables. An improvement of the 
bound on the /1 — > e + 7 decay rate by a factor of ~ 10 will make the observation of pseudo-Dirac neutrinos 
at the LHC completely impossible and the effects on the electroweak precision observables negligibly small, 
unless the RH neutrinos have additional (flavour-conserving) couplings to the Standard Model particles, as 
like in theories with extra U(l) local gauge symmetry under which the pseudo-Dirac neutrinos are charged, 
or the TeV scale type III see-saw scenario. 



V. PREDICTIONS FOR THE (/3/3)o„-DECAY 



In the scenario we are considering, the (/3/3)o„-decay effective Majorana mass |<m>|, which controls the 
(/?/3)oi/-decay rate, receives a contribution from the exchange of the heavy Majorana neutrino fields N k |12) . 
which may be not negligible for "large" neutrino Yukawa couplings. One has, in general: 



<m> 3* 



J2 U 2 el m,-^ F(A, M k ) (RV) 2 ek M k 



k=l 



(5.1) 



where 7 the function F(A, M k ) depends on the N k masses and the type of decaying nucleus (A, Z). For M k = 
(100— 1000) GeV, one can use the rather accurate approximate expression for F(A,M k ) (T3]: F{A,M k ) = 
(M a /M k ) 2 f(A), where M a w 0.9 GeV and f(A) depends on the decaying isotope considered. For, e.g., 76 Ge, 
82 Se, 130 Te and 136 Xe, the function f(A) takes the values f(A) ^ 0.079, 0.073, 0.085 and 0.068, respectively. 
In the case of 48 Ca, f(A) has a smaller value [13] : /( 48 Ca) = 0.033. Using eq. (2.16) we can write the heavy 
Majorana neutrino exchange contribution to |<m>| in a simplified form: 



<m> 1N = - 



This contribution, as we will show below, can be even as large as |<m> Ar | 
(IH) light neutrino mass spectrum. 8 



(5.2) 

0.2 (0.3) eV in the case of NH 



7 Let us note that the interference between the contributions due to the light and heavy Majorana neutrino exchanges is not 
suppressed because both contributions are generated by the weak interaction involving currents of the same (V — A) structure. 
The interference term of interest would be strongly suppressed if the heavy Majorana neutrino exchange is generated by (V+A) 
currents |1 21 . 

8 Note that in the approximation we use for (RV) e k one has ^2 k —i (RV)^k Mk = 0. The heavy Majorana neutrino exchange 



contribution <m> to |<m>| is not zero due to the nontrivial dependence on M k of the function F(A, M k ), see eq. (5.1 1 
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In Fig. [3] we show the ratio between the total effective Majorana mass |<m>| given by eq. (5.1) and 



tdl _ 



the "standard" contribution due to the light Majorana neutrino exchange (see, e.g., [36], EZ]) : |<m> s 
I Xh=i m i\- I n this plot we considered a nuclear matrix element factor f(A) corresponding to 76 Ge, 
although the conclusions are analogous for other nuclei. The scan of the parameter space was done in the 
same way as in Fig. [T] selecting just the points which are in agreement with the present experimental bound 
on B{ji — > e + 7). In Fig. [4] we show the range of the possible values of |<m>| as function of |(i?V) e i| for 
Mi = 100 GeV and z < 1CT 2 , in the case of (/?/3) „-decay of 76 Ge. 

The results of the analysis illustrated graphically in Figs. [3] and [4] demonstrate that RH neutrinos can 
significantly enhance the rate of (/3/3)o;/-decay, without this being in conflict with the present upper bound 
on the fi — > e + 7 decay rate. For example, |<m>| can be as large as 0.2 (0.3) eV, if z = 10~ 3 (10 -2 ) and 
Mi 1=3 100(1000) GeV. In contrast, in the limit z — > 0, the RH neutrinos behave as a Dirac pair and hence 
do not contribute to the (/3/3)oi;-decay amplitude to leading order. 

Let us recall that in the scheme with two heavy Majorana neutrinos we are considering, the lightest neutrino 
mass is zero. This implies that for the NH and IH light neutrino mass spectra we have for the "standard" 
contribution to |<m>|, respectively (see, e.g., [37]): |<m> std | & 0.005 eV and 0.01 eV & |<m> std | & 0.05 
cV. Thus, if it is established that the light neutrino mass spectrum is hierarchical, a value of | < m > | > 0.05 eV 
would signal the presence of a contribution to |<m>| beyond the "standard" one. In the scheme considered, 
|<m>| can be by a factor up to ~ 100 (~ 10) larger than the maximal value of |<m> std | predicted in the 
case of NH (IH) light neutrino mass spectrum. 

It should be noted also that the predicted value of | < m >| in the cases of the NH (IH) spectrum exhibits 
a strong dependence on the PMNS parameters and especially on the Majorana phase |31j CX21 — «3i (o^i). 
If we have |<m> std | = |<m> N |, i.e., if the "standard" contribution to |<m>| is of the same order as the 



contribution due to the exchange of the heavy Majorana neutrinos Ni l2 (see eq. (5.2)), |<m>| will depend 



also on the phase u) (see eq. (2.15)). 

Finally, if for a given decaying nucleus ( A\ , Z\ ) and certain values of the parameters of the problem there 
is a strong mutual compensation between the two contributions <m> std and <m> N in |<m>| and we 
have |<m>| <C |<m> std |, |<m> N |, similar cancellation will not happen, in general, for another decaying 
nucleus (A2,Z 2 ) due to the dependence of <m> N on (A, Z) [T2]. In the scheme considered, in which only 
hierarchical light neutrino mass spectrum is possible, and in view of the planned sensitivity of the next 
generation of (/3/3)o„-decay experiments, this observation has practical importance if the light neutrino mass 
spectrum is with inverted hierarchy. In a more general context in which the 3rd heavy Majorana neutrino is 
relevant for the see-saw mechanism and the light neutrino mass spectrum can be with partial hierarchy or 
even of the quasi-degenerate type (see, e.g., [551 [35]), the observation made above regards both types of light 
neutrino mass spectrum - with normal ordering and with inverted ordering. Indeed, <to> n in the general 
case can be written in the form: 

3 M 2 

< m > N < A - -f(A) F N , F N = (RV) 2 ek ^ , (5.3) 

k=l k 

where F N does not depend on A. Suppose that for a given isotope (A 1} Z\) the two terms in |<m>|, <m> std 
and <m> N ' Al , mutually compensate each other so that \<m> Al \ <C |<m> std |, |<m> N ' Al |. That would 
imply that 

3 

<m> N ' A ^f(A 1 )F N ^U>„ (5.4) 

i=l 

In this case the effective Majorana mass corresponding to an isotope (A2, Z2) will be given by: 

KM) 



\<m> A2 \ = 



1 - 



(5.5) 



If, e.g., the cancellation between the two terms in <m> Al occurs for 48 Ca (for which /( 48 Ca) = 0.033), 
it will not take place for e.g., 76 Ge, 82 Se, 130 Te and 136 Xe since for these isotopes f{A 2 ) = 0.079, 0.073, 
0.085 and 0.068, respectively. Actually, for 76 Ge, 82 Se and 130 Te the factor |1 - f{A 2 )/f(A 1 )\ = 1.39, 1.21 
and 1.58, so we will have a somewhat larger |<m>" 42 | than the "standard" one |<m> std |, while for 136 Xe 
the indicated factor is 1.06 and thus |<m>" 42 | = |<m> std |. If, however, the cancellation between the two 
terms in <m> Al takes place for, e.g., one of the nuclei 76 Ge, 82 Se, 130 Te and 136 Xe, for which the function 
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f{A) has rather similar values, |<?7i> A2 | for the other nuclei will be suppressed with respect to |<m> std | 
to various degrees. For instance, if the cancellation is operative for 136 Xe, |<m> j42 | for 76 Ge, 82 Se, 130 Tc 
and 48 Ca will be suppressed with respect to the "standard" contribution |<7n> std | by the factors 0.16, 0.07, 
0.25 and 0.51, respectively. 

In the case with two heavy Majorana neutrinos we are considering and for IH light neutrino mass spectrum, 
the condition for an exact cancellation between <m> std and <m> N ' A can be easily derived in terms of the 



basic parameters of the scheme. Using eqs. (2.18), (2.19), (5.1) and (5.2) we can write |<m>| as: 
|<m>| ~ \ mi U^{l~K) + m 2 U 2 e2 {l + K) + 2i^n^rnr 2 K (U e2 U el ) I , 



where K is given by: 



K 



9 9 

z y v 



-f(A)e 



(5.6) 



(5.7) 



If we require |<m>| ~ 0, the factor K, which depends on the see-saw parameters y, z, uj and Mi, is 
expressed only in terms of the neutrino oscillation parameters: 



K 



cos 20i2 + i sin 20i 2 cos ^ 
1 + sin 20i2 sin ^ 



If the Majorana phase «2i takes a CP-conserving value we get: 

COS 2012 



K 



for a 2 i = (2k + l)?r , fc = 0,±l, 



1 + r] k sin (20i2 

e ^V k 6l2 for 

where rjk = (— l) fe - In the case of <m> Al = we obtain for |<m> A 



K ~ e 2ink912 for a 2 i=2fc7r, fc = 0,±l,..., 



\<m>~ 2 \ = 



1 - 



f(A 2 ) 



f{A 1 ) 



13 



1 - sin 20i2 siu 



2 a 2l\ 2 



(5.8) 

(5.9) 
(5.10) 

(5.11) 



The condition (5.8) (or (5.9) and (5.10)) for <m> A = strongly constrains the phase uj and the size 
of K. In order to be satisfied, the condition of cancellation between <m> std and 
correlat ion betw een the values of the see-saw parameters y, z, the phase uj and Mi 



<m> N,A requires a 
and the neutrino 



m 2 , the solar neutrino mixing angle 0i2 and the Majorana phase a2i- A priori, such 



a correlation seems highly unlikely, making the cancellation in the effective Majorana mass in the scheme 
considered appear unprobable. Thus, if the light neutrino mass spectrum is with inverted hierarchy, within 
the see-saw scheme considered with two heavy Majorana neutrinos, it appears quite unlikely that if, for 
example, the GERDA III experiment with 76 Ge observes a positive (/3/3)oi/-decay signal, another (/3/3)o v - 
decay experiment which uses a nucleus different from 76 Ge, will not see a signal due to a strong suppression 
of the effective Majorana mass caused by the cancellation under discussion for that nucleus, and vice-versa. 

The dependence of the interplay between <m> std and <m> N 
due to the exchange of light and heavy Majorana neutrinos, on 
(IH spectrum) or 0:21 — «3i (NH spectrum) is illustrated in Fig.j5|for Mi = 100 GeV and y = 0.01. The solar 



i.e. between the contributions to |<m>| 
z, uj and the effective Majorana phase a 2 \ 



and atmospheric neutrino oscillation parameters, 



Ar 



and 



7 23> 



), respectively, are fixed to 



their corresponding best fit values, and we have set 0i3 = 0.2 and S — 0. Notice that the plots showing the 
correlation of |<m>| and the Majorana phase are symmetric with respect to a2i or (021 — 0:31) equal to 
7r/2 and 37r/2 if the phase u> takes the values uj — 0,7r/2,7r. As Fig. [5]shows, we can have |<m>| 0.01 eV 
in the case of NH spectrum for a relatively large range of values the Majorana phase 021 — 0131 and certain 
values of the phase uj: for, e.g., uj — 7r/4, we get |<m>| ^ 0.01 eV for 5 ^ 0:21 — a^i ^ Air if z — 10 -3 . 
Actually, for the indicated values of z and a 2 i — 031 we have |<m>| ^ 0.01 eV for any value of uj from the 
interval [0, 2ir]. The interplay between the Majorana phase and uj can induce also a minimum of |<m>| for 
a certain value of the degeneracy parameter z. For example, in the case of the IH spectrum, for a 2 i = 
and uj — 7r/3, the predicted effective Majorana mass can be smaller than 0.03 eV: we have for the different 
nuclei considered 0.007 eV < |<m>| < 0.03 eV if 2 x 10~ 4 < z < 10" 3 . 
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FIG. 3: The ratio between the effective Majorana mass [<m>| and the "standard" contribution |<m> std | for 76 Ge 
in the cases of NH (left panel) and IH (right panel) light neutrino mass spectrum, for Mi — 100 GeV and i) y = 0.001 
(blue o), ii) y — 0.01 (green +). 
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FIG. 4: The effective Majorana mass \<m>\ vs \ (RV) e i \ for 76 Ge in the cases of NH (left panel) and IH (right panel) 
light neutrino mass spectrum, for Mi = 100 GeV and i) y = 0.001 (blue o), ii) y — 0.01 (green +). The gray markers 
correspond to |<m> stcl |. 



VI. THE (/3/3V-DECAY EFFECTIVE MAJORANA MASS AND B(fi e + 7) 



We now combine the information on the see-saw parameter space that we obtained from the lepton flavour 
(jit — > e + 7) and lepton number ((/3/3)oy) violating processes studied in the previous sections. In the simple 
extension of the Standard Model considered so far, i.e. with the addition of two heavy RH neutrinos iVi 
and N 2 at the TeV scale, which behave as a pseudo-Dirac particle and at the same time are responsible 
for the generation of neutrino masses via the see-saw mechanism, a sizable (dominant) contribution of Ni 
and N 2 to the (/3^)o J/ -decay rate would imply a "la rge" effect in the mu on ra diati ve decay rate. Indeed, if 



|<m>| = |<m> N |, where <m> N is given in eq. (|5.2[), using eqs. (|3.1[) and (|3.5[) it is easy to show that, 



given the splitting 10 < z <§; 1 between Mi and M2, |<m>| 
—> e + 7 decay branching ratio. More explicitly, we have: 



\<m> I can be directly related to the 



B(p -4- e + 7) = 



3a r 



64tt 



\G(0)~G(X)\ 2 \r\ 



Ml 



<m> 



N|2 



(6.1) 
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FIG. 5: Upper panels: the dependence of the effective Majorana mass |<m>| on the Majorana phase Q21 — Q31 
(left side), 0:21 (right side) for M\ = 100 GeV, to = tt/4, y — 0.01 and z = 10 . The gray dash-dotted line shows 
j<ra> std j. Lower panels: |<m>| versus the degeneracy parameter z for Mi = 100 GeV, uj = 0, y = 0.01 and 
ct2\ — 0. The left (right) panels correspond to a NH (IH) light neutrino mass spectrum. 



where r = (Uu2, ~~ *\/ TO 3/ TO 2 U^)/(U e 2 — i^/m^/m2 U e ^) for the NH mass spectrum. As was pointed out 
earlier, the corresponding expression for the case of IH spectrum is obtained by replacing TO2,3 - ► Wi,2 
and U a 2, a 3 —> U a \,ai [p. — e,fi). For the NH (IH) neutrino mass spectrum one has: 0.5 < \r\ < 30 
(0.01 < \r\ < 5). Therefore, for Mi « 100 GeV and NH spectrum, the predicted rate of the /i — ?► e + 7 decay 
can even be larger by up to one order of magnitude than the upper bound on B(fi — » e + 7) if a positive 
signal in the current experiments searc hing for (/3/3)oi/-decay is detected, implying |<m>| <~ 0.1 eV. 

A lower bound on B(/x — > e + 7) in (6.1 ) can be derived for both types of light neutrino mass spectrum. 
In the case of NH spectrum, such lower bound is within the sensitivity of the MEG experiment, provided a 
(/3/?) iy-decay corresponding to |<m>| ^ 2 x 10~ 2 eV is observed. 

The analytic relation between B(fi — ¥ e + 7) and |<m>| in eq. (6.1) is confirmed by the results of 
numerical computation, reported in Fig. [6} The plot shows the correlation between the fi — > e + 7 branching 
ratio and the effective Majorana mass in the case of "large" couplings between the RH (pseudo-Dirac pair) 
neutrinos and charged leptons in the Lagrangian (2.3). The effective Majorana mass |<m>| was computed 
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for z = 10 using the general expression (5.1). The (/3/3)oiy-decay nucleus was assumed to be Ge. The 
neutrino oscillation parameters are taken, again, within the corresponding 3er experimental intervals reported 
in Table 1. The Majorana phase 021 (031 — 021) and the phase uj in the IH (NH) case were varied in the 
intervals [0, 47r] and [0, 2n], respectively. The neutrino Yukawa cou plin g takes values y < 0.1. The correlation 
between B(/i — > e + 7) and |<m>| = |<m> N | reported in eq. (6.1 1 is satisfied for values y > 0.01. This 
is in agreement with Figs. [3]and|4j where it is shown that a signal compatible with the GERDA sensitivity 
reach is possible, provided y > 10~ 3 , for both types of neutrino mass spectrum. Moreover, in the case of IH 
light neutrino mass spectrum, such correlation depends strongly on the value of the Majorana phase 0121. 
Indeed, for Mi = 100 (1000) GeV and y = 0.01 (0.1) we expect that the MEG experiment (34] is able to 
measure the fi — > e + 7 decay rate (see Fig. [2]). If lepton flavour violation is discovered by MEG, according 
to eqs. (5.2) and ( |3.10| , a positive signal detected by GERDA II, i.e. |<m>| = \<m> N \ > 0.1 cV, implies: 
IO- 3 (lCT 7 ) < z(l + 0.94sin(a 2 i/2)) < 4 x 10" 3 (4 x 10" 2 ). In the case of M x = 100 GeV and z = 10~ 3 , 
used to obtain Fig.|6j we would expect, in general, positive signals to be observed in both MEG and GERDA 
II experiments if a 2 \ = 0,7r; in the case of 021 = 37r, the (/3/?)oj, and /1 — >• e + 7 decays are predicted to 
proceed with rates below the sensitivity of these two experiments. 

We note, however, that it is not possible to get independent constraints on the degeneracy parameter z 
and the Majorana phase from the data on (P(3)o v and [i — > e + 7 decays. Finally, we notice also that the 
strong correlation exhibited in Fig. [6] is a consequence of the constraints imposed by the neutrino oscillation 
data on the type I see-saw parameter space in the case investigated by us. 



VII. CONCLUSIONS 



We have analyzed the low energy implications of a type I see-saw scenario with right-handed (RH) neutrino 
masses at the electroweak scale and sizable charged and neutral current weak interactions. This class 
of scenarios have the attractive feature that the RH neutrinos could be directly produced at the Large 
Hadron Collider, thus allowing to test in collider experiments the mechanism of neutrino mass generation. 
Furthermore, and in contrast to the high-scale see-saw mechanism, the rates for the rare leptonic decays are 
unsuppressed in this scenario, which opens up the possibility of detecting signatures of new physics with 
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experiments at the intensity frontier. 

Present low energy data set very stringent constraints on this scenario. Namely, reproducing the small 
neutrino masses requires, barring cancellations, that two of the heavy RH neutrinos must form a pseudo- 
Dirac pair even in the case when there is no conserved lepton charge in the limit of zero splitting at tree 
level between the masses of the pair. Besides, reproducing the experimentally determined values of the 
low energy neutrino oscillation parameters (mixing angles and neutrino mass squared differences) fixes the 
weak charged current (CC) and neutral current (NC) couplings of t he h eavy Majorana neutrinos to the 
charged leptons and W ± and light neutrinos and Z°, (RV)ek (see eqs. (2.3) and (2.4)), up to an overall scale 
which can be related to the largest eigenvalue y of the matrix of neutrino Yukawa couplings. This allows 
to derive explicit expressions for the rates of the lepton flavour violating (LFV) charged lepton radiative 
decays, fi — > e + 7, t — > e + 7, t — > /i + 7, in terms of the low energy (in principle measurable) neutrino 
mixing parameters (including the Dirac and Majorana CP violating phases), the neutrino Yukawa coupling 
y and the RH neutrino mass scale. Using the present constraint on the rate of the process /x — ¥ e + 7 we 
have obtained an upper bound on the coupling y under the assumption that the RH neutrino mass scale is in 
the range (100 -j- 1000) GcV. Our analysis shows that the restrictions on this scenario from the data on the 
neutrino mixing parameters and the upper bound on the /i — > e + 7 decay rate imply that the CC and NC 
couplings |(i?U)tt| of the heavy RH Majorana neutrinos are too small to allow their production at the LHC 
with an observable rate. Other lepton flavour violating processes, such as /i — e conversion in nuclei, give 
complementary constraints on the parameter space of this model, which will be discussed elsewhere [35]. 

We have also analyzed the enhancement of the rate of neutrinoless double beta ((/3/3)oi/-) decay induced 
by the RH neutrinos. We have shown that even after imposing the restrictions on the parameter space 
implied by the data on the neutrino oscillation parameters and on the LFV charged lepton radiative decays, 
the contribution due to the exchange of the RH neutrinos in the (/3/3)oi,-decay amplitude can substantially 
enhance the (/3/3)o„-decay rate. As a consequence, the latter can be in the range of sensitivity of the GERDA 
experiment even when the light neutrinos possess a normal hierarchical mass spectrum. Finally, the rate 
of the \i — > e + 7 decay, generated by the exchange of the heavy RH neutrinos can naturally be within 
the sensitivity of the MEG experiment. Thus, the observation of (/^o^-decay in the next generation of 
experiments which are under preparation at present, and of the fi — > e + 7 decay in the MEG experiment, 
could be the first indirect evidence for the TeV scale type I see-saw mechanism of neutrino mass generation. 
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